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ABSTRACT 

We report the discovery of three pulsars whose large dispersion measures and 
angular proximity to Sgr A* indicate the existence of a Galactic center popula- 
tion of neutron stars. The relatively long periods (0.98 to 1.48 s) most likely 
reflect strong selection against short-period pulsars from radio-wave scattering 
at the observation frequency of 2 GHz used in our survey with the Green Bank 
Telescope. One object (PSR J1746-2850I) has a characteristic spindown age of 
only 13 kyr along with a high surface magnetic field ~ 4 x 10 13 G. It and a 
second object found in the same telescope pointing, PSR J1746-2850II (which 
has the highest known dispersion measure among pulsars), may have originated 
from recent star formation in the Arches or Quintuplet clusters given their an- 
gular locations. Along with a third object, PS R J1745-2910, a nd tw o similar 



high-dispersion, long-period pulsars reported by iJohnston et al.l (120061 ). the five 
objects found so far are 10 to 15 arc min from Sgr A*, consistent with there being 
a large pulsar population in the Galactic center, most of whose members are un- 
detectable in relatively low-frequency surveys because of pulse broadening from 
the same scattering volume that angularly broadens Sgr A* and OH/IR masers. 

Subject headings: pulsars: general — pulsars: individual (J1746-2850I, J1746- 
2850II, J1745-2910) 



1. Introduction 



The scientific pay-offs from finding pulsars orbiting near Sgr A* are potentially very 
high and fall into three main categories. Measurements of relativistic effects through timing 
of pulsars in tight orbits around Sgr A* would provide methods for better constraining the 



mass of the central bla c k hole and even esti mating; its spin (e.g. lLaguna fc Wolszczanl 119971 . 
Wex fc Kopeikinl Il999i . iPfahl fc Loebl 120041 ) . Timing measurements will also characterize 
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the distribution of dark matter near Sgr A* either in the form of a cluster of black holes 
and neutron stars or in a smoothly distribute d volume containing dark-matter particles 
( IBertone &: Merrittl 120051 . IWeinberg et al.ll2005l ). The spatial, age and period distributions 
of pulsars near Sgr A* will help describe the stellar population in the region and discriminate 
between hypotheses attempting to explain the presence of the central cluster of young massive 
stars: stellar collisions and mergers, migration, and a past episode of intensive star formation. 
Finally, using pulsars to probe the scattering region around Sgr A* will lead to refinement 
of electron density models for the inner Galaxy. 

There are at least three stellar clu sters in the Galact ic center region containing massive 
stars, the progenitors of radio pulsars. iGhez et al.l (120031 ) present evidence that the central 
cluster around Sgr A* includes early- type stars. Since stars with masses of about 8 — 20 M Q 
are neutron star (NS) progenitors, it is plausible that a considerable NS population exists 
as well. This leads to the con clusion that a sizable fraction of the NSs would be active radio 
pulsars. iPfahl fc Loebl (12004 ) estimate tha t there are ~ 10 — 10 00 active pulsars orbiting 



Sgr A* with periods of less than 100 years. ICordes fc Lazio Jl997h show that some of these 
pulsars will be detectable at the distance of the Galactic center, even at the high frequencies 
needed to mitigate pulse broadening from scattering. 

Apart from the central star cluster around Sgr A*, the Arches and Quintuplet clusters 
are the densest stellar associations within 0.5° of the Galactic center. Both clusters are 
composed mainly of young massive stars whose supern ova explosions will have produced 
neutron stars young enough to be active radio pulsars. IStolte et al.l (120081 ) show that the 
radial distance from the Galactic center to the Arches cluster is 62 ± 23 pc if the cluster is 
in a circular orbit around Sgr A*, but its distance may be as large as 200 pc. 

Of the nearly 2000 pulsars known up to now, only two are within 1° of Sgr A*: PSR J1746— 2856 
and PSR J1745— 2912, which were discovered by I Johnston et al.l (120061 ) in a 3.1 GHz survey 
with the Par kes telescope and hav e dispersion measures (DMs) of 1168 and 1130 pc cm -3 , 
respectively. iJohnston et al.l (120061 ) also report the absence of detections by a 8.4 GHz Parkes 
survey, which indicates there are not many pulsars with flat enough spectra to be detected at 
high frequencies. The non-detections are also consistent with electron- density models that 
predict that scatteri ng at frequencies up t o ~ 10 GHz strongly affects the detectability of 
pulsars near Sgr A* (ICordes fc Laziolll997l ). 



As part of a larger program to probe the Galactic center pulsar population and ionized 
gas environment, we have surveyed the inner 0.5° around Sgr A* at 2 GHz with the Green 
Bank telescope. Pulse broadening due to scattering varies with frequency approximately as 
t s oc z/~ 4 , so the broadening expected at 2 GHz is ~ 5 times larger than at 3.1 GHz. This 
effect is partially mitigated by the fact that most pulsars have moderate to steep spectra 
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and are more easily detected at lower frequencies when scattering is not a factor. 



2. Survey Parameters and Data Processing 

The survey was carried out in September 2007 with the Green Bank telescope, using a 
center frequency of 1.95 GHz and the SPIGOT backend with a bandwidth of 800 MHz. The 
lower 200 MHz of the receiver band was excluded from the analysis due to the presence of 
severe radio frequency interference. The effective analyzed bandwidth was 600 MHz divided 
into 768 channels. Data were taken with a sampling time of 82 fis in two polarizations, 
which were summed before further processing. The grid of survey pointings around Sgr A* 
is shown in Fig. [1] overlaid on a 330 MHz image of the Galactic center region. We observed 37 
grid pointings for one hour each, but data from 8 pointings were not usable due to technical 
issues and these pointings are omitted from the figure. 

Data processing was performed at the Cornell Center for Advanced Computing (CAC) 
using the Cornell pulsar search code3 and the Presto packag^E The full-resolution data were 
searched with 1245 evenly spaced trial DM values in the range 10 - 2000 pc cm -3 , and 155 
evenly spaced values in the range 2000 - 2500 pc cm -3 . Data were also searched over 1000 
values from 100 to 2100 pc cm -3 with time resolution of 0.65 ms and 300 values from 2100 
to 3000 pc cm -3 with 1.3 ms resolution. A Fast Fourier Transform was performed on the 
resulting time series and the power spectrum was searched for periodic signals by summing 
up to 16 harmonics and applying a harmonic sum threshold of 6a. 

Time series were also searched for individual dispersed pulses using two algorithms, 
one that applies simple templates (Cordes & McLaughlin 2003) and another that searches 
for clusters of related points (Cordes et al. 200 4) using a friends-o f-friends group-finding 



algorithm (see Fig.l in lHuchra fc Gellerlll982l ; see lDeneva et al.ll2008l for detailed description 



of application to single pulse searching). The matched filtering approach smoothes each 
time series with rectangle functions of widths 2 n samples, where n = — 7, and searches for 
events above a 5a threshold. It is most sensitive to pulses whose widths are close to one of 
the boxcar filter widths, for a maximum width of 10 ms. The cluster algorithm identifies 
individual samples above a 3a threshold and then combines contiguous points into events to 
which the larger 5a threshold is applied. 



^ttp : / / arecibo . tc . Cornell . edu/PALFA/ 
2 http : //www . cv . nrao . edu/^sransom/presto 
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Fig. 1. — The grid of 2 GHz survey po intings from 2007 overlaid on a 0.33 GHz image of 
the Galactic center (lLaRosa et al.ll2000i ). Circles correspond to the 5.8' FWHM beam size 
of the Green Bank telescope at 2 GHz. Diamonds denote the positions of J1746— 28501, 
J1746-2850II, and J1745-2910, and triangles denote the positions of J1746-2856 and 
J1745-2912. Crosses show actual position uncertainties. For J1746-2850I, J1746-2850II, 
and J1746— 2856 the position uncertainties are smaller than the marker size. Squares show 
the positions of the Arches and Quintuplet clusters. 
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3. Survey Results 

Table [I] summarizes the properties of the three discovered pulsars. PSR J1746— 28501 
and PSR J1746— 2850II were discovered in the same survey pointing and were confirmed 
with the Green Bank Telescope in June 2008; subsequent monthly timing observations at 
2 GHz are ongoing. We observed both pulsars at 1.5 and 4.8 GHz with the dual goals 
of estimating their spectral indices and improving their position measurements. We also 
observed PSR J1946— 28501 at 9 GHz in order to confirm that its relatively flat spectrum 
extends to higher frequencies. PSR J1745— 2910 was discovered in February 2009 and will 
be timed starting in June 2009. Profiles of the three pulsars are shown in Figure [2J 

We estimated the period-averaged flux density for each pulsar by scaling from signal- 
to- noise ratio of the folded pulse profile {S/N) { using the calculated radiometer noise level 
for the sum of two polarization channels: 

s _ (S/N) pio{ T sys 



G-s/N poY Av T obs /n b in' 

and we have used T sys of 32 K, 27 K, 19 K, and 27 K at 1.4, 2, 4.8, and 9 GHz, including sky 
background contribut ions in the Galactic center direction of 12 K at 1.4 GHz scaled as v~ 2 ' 5 
faeich fc Reichlll988h . The gains G are 2.0 K Jy" 1 , 1.9 K Jy" 1 , 2.0 K Jy" 1 , 1.8 K Jy" 1 , 



respectively at the four frequencies. The number of bins in the pulse profile nbi ns — 128. The 
total bandwidth Av = 600 MHz at 1.4 and 2 GHz, and 800 MHz at 4.8 and 9 GHz. 

No isolated dispersed pulses were detected above a 5a threshold in any of the 2 GHz 
survey pointings. For a 100 ms pulse width (comparable to the observed pulse widths from 
Tabled]), we derive an upper limit on the flux density of single pulses 

S max = - M mTsy ; « 6.5 mJy (2) 
GyiVpoi Av w 



3.1. PSR J1746-2850I 

A timing solution for J1746-2850I was obtained from observations with the GBT made 
at 24 epochs between June 2008 and August 2009 and fitted for period P, period derivative 
P, DM, and sky coordinates. The long period and large period derivative indicate that 
this is a young object with a characteristic age r c = P/2P = 13 kyr and large surface 
magnetic field, 3.8 x 10 13 G. Assuming a power-law spectrum S v oc u a , we found that J1746- 
28501 has a relatively flat spectrum, with a = —0.3. The mean spectral index for normal 
pulsars is —1.5 and fewer than 10% have a > —0.5 (Lorimer et al. 1995). Of the two 
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Table 1: Three New Pulsars Toward the Galactic Center 



Parameter 


J1746-2850I 


J1746-2850II 


J 1745-29 10 


RA (J2000) 


yjh 46 m 06.6 s (2) 


yjh 46 m Q3.7 S (1) 


yjh 45m 16^(34) 


DEC (J2000) 


-28° 50' 42"(5) 


-28° 49' 19"(21) 


-29° 10'(3) 


P(s) 


1.0771014910(4) 


1.478480373(2) 


0.982 


P (s/s) 


1.34311(2) xl0~ 12 


1.27(6) xl0~ 14 




DM (pc cm- 3 ) 


962.7(7) 


1456(3) 


1088 


Age (Myr) 


0.013 


2 




B(G) 


3.8 x 10 13 


2.8 x 10 12 




E (erg/s) a 


4.24 x 10 34 


1.47 x 10 32 




W 8 .g (50%,ms) fe 


50 






W4.8 (50%,ms) 


30 


31 




W 1 . 95 (50%,ms) 


45 


130 


54 


W 1A (50%,ms) 


100 


145 




Sg.g (mJy) 


0.4 






S 4 . % (mJy) 


0.5 


0.1 




51.95 (mJy) 


0.6 


0.2 


0.2 


51.5 (mJy) 


0.8 


0.4 




Spectral index a 


-0.3 


-1.1 




(S u oc u a ) 









a Assuming a 1.4 M Q neutron star with a 10 km radius and moment of inertia I = 10 45 g cm 3 . 
b Subscripts refer to observing frequencies in GHz. 
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currently known r adio-emitting magn etars, XTE1810— 197 has a = —0.5 in the frequency 

l|2006h . t or rising spectrum 



range 0.7-42 GHz (Camilo et al 



over 1.4 — 6.6 GHz ( I Camilo et al 



2003). The flat spectrum of PSR J1746-2850I in addition 



to its spin parameters suggest that this object belongs to a class of young neutron stars 
bridging the gap between magnetars and canonical radio pulsars. 



3.2. PSR J1746-2850II 

J1746— 2850II was timed simultaneously with J1746— 28501 because the two pulsars are 
contained within the same beam of the GBT at 1.4 and 2 GHz. It has the highest dispersion 
measure of any pulsar known to date, indicating that of the five pulsars currently known 
within 1° of Sgr A*, J1746-2850II is likely closest to Sgr A* in radial distance. It is an old 
pulsar, with a characteristic age r c = 2 Myr. While most of the timing observations were 
made at 2 GHz, the receiver was unavailable in July and August 2008 so observations were 
made at 1.5 GHz. J1746— 2850II is significantly scattered at this frequency (Fig. [2]). While 
it was detected during a 2.25 h observation in June 2008, it was too weak to detect in our 
routine 1-hr observations at 1.5 GHz. Consequently, its timing solution has a larger position 
uncertainty than that of J1746-2850I. 

The pulse profile of J1746— 2850II at 4.8 GHz is symmetric and Gaussian-like without 
a discernible scattering tail (Fig. [2]). We simulate the effect of scattering broadening by 
convolving a Gaussian matching the 4.8 GHz profile with an exponential depending on the 
scattering time, t sc . We compare the results to the 1.5 and 2 GHz profiles using least-squares 
fitting and obtain t sc pa 140 ms at 2 GHz and 266 ms at 1.5 GHz. 



3.3. PSR J1745-2910 

PSR J1745-2910 has a period of 0.982 s and DM of 1088 pc cm -3 . Its nominal position 
is 12' from Sgr A* and ~ 8' from J1745— 2912. Due to the large uncertainties in the positions 
of both pulsars, they may be between 3' — 14' apart. Using Eqn. [TJ we estimate the flux 
density of J1745-2910 to be ~ 0.2 mJy at 2 GHz. 



3.4. PSR J1746-2856 and PSR J1745-2912 



The two pulsars discovered at 3.1 GHz by iJohnston et al.l (120061 ) are in our 2 GHz 
search area (Fig. [T]). PSR J1746— 2856 was blindly detected at up to ~ 6' from its position. 
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PSR J1745-2912 was not detected either blindly or by folding data from search pointings 
near its position with its period of 187 ms. The reported scattering time for this pulsar is 
25 ± 3 ms at 3.1 GHz. If we assume a Kolmogorov scattering spectrum with a dependence 
of the scattering time on frequency of r sc oc /~ 40 , the scattering time for PSR J1745-2912 
at 2 GHz is 144 ms. This is close to its period and therefore we attribute the non-detection 
to scattering broadening. 



4. The Galactic Center Environment 

Scattering from the dense regi on near Sgr A* has bee n modeled as both a thin screen and 



as an extended scattering volume. lLazio fc Corded (119981 ) estimate that the scattering screen 



aroun d the Galactic center is 133lgo° P c fr° m Sgr A*. The NE2001 model (ICordes fc Lazio 



20031 ) uses an ellipsoidal scattering volume with a scale height of 26 pc and characteristic 
radius of 145 pc that is slightly offset from Sgr A*. The expected DM of a pulsar at the 
location of Sgr A* is ~ 1600 pc cm -3 with a corresponding pulse-broadening time ~ 400 s 
or ~ 2000 s at 1 GHz for the screen or extended models, respectively. Pulsars on the near 
side of Sgr A* will have smaller values of DM and much smaller pulse broadening times. 

All five pulsars within 1° of Sgr A* have high DMs that are too large to be accounted 
for by the Galactic disk components of the NE2001 electron density model. If the pulsars 
are part of a disk population in the foreground of Sgr A* the large DMs might be due 
to an unmodeled e xcess of interve n ing io nized gas associated with foreground HE regions. 



However, the two iJohnston et al.l (120061 ) pulsars are far in angular separation from one 
another and from the three pulsars we have found, making unlikely the possibility of a single 
HII region contributing to their large DMs. Similarly unlikely is the coincidence of five out 
of five pulsars being affected by separate HII regions. 

In the absence of any unmodeled HII regions, the NE2001 model places all five pulsars 
within ~ 100 pc of Sgr A* using the DM values (Fig. [3]). However the predicted pulse 
broadening times for the pulsars are several orders of magnitude larger than the measured 
value for J1746— 2850II and the upper bounds on the other objects. This signifies that 
the near-side boundary and/or density of the idealized Galactic-center compo nent of the 



model are incorrect or that the scattering volume is patchy, as suggested by lLazio et al. 



(119991 ). In either case, like any distance calculated using DM, the distances of the five pulsars 
from Sgr A* are highly model dependent and cannot be refined to better than ~ ±100 pc. 
However, unless there is an electron density component in the foreground disk (or spiral arm) 
of the Galaxy, the large DMs imply that the pulsars are no more than about 200 pc from 
Sgr A*. Galactic center pulsars will be key for improving the Galactic center component of 
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the next electron density model, NE2008 (Cordes et a!., in preparation). 



4.1. The Arches and Quintuplet Clusters 



Could J1746— 28501 have been born in either the Arches or Quintuplet clu ster? Both 



cluste rs contain numerous massive stars, which are neutron star progenitors. iKim et al. 



(120061) find the present day mass function of the Arches cluster to extend to at least 40 M Q , 
while iFiger et al.l (120021 ) find the Arches cluster to contain stars exceeding 100 M with the 
age of the cluster estimated to be 2 ± 1 Myr, old enough for neutron stars to have formed. 
Current estimates for the maximum progenitor mass that will produce a neutron star are 
20 M for solitary progenitors a nd as high as 50 — 80 M ffi for binary progenitors due to 
mass loss in Roche lobe ove rflow (IBelczynski fc Taamll2008l ). The Quintuplet cluster has an 
estimated age of ~ 4 Myr (IFiger et al.lll999l ) but contains the Pistol Star, whose mass of 
150 M points to an age no larger than 2 Myr. These estimates indicate that the Quintuplet 
cluster is also old enough to have already formed neutron stars. 

The current position estimate for J1746— 28501 places it within 2' of the Quintuplet 
cluster (Fig. [I]). If the cluster and the pulsar are roughly 8.5 kpc from Earth, this angle 
corresponds to a distance of 5 pc. A transverse velocity of 375 km/s would have allowed the 
pulsar to reach its present location from the Quintuplet cluster within 13 kyr. Our double 
discovery in a single pointing suggests that the Arches and Quintuplet region may have a 
relative excess of pulsars. 



4.2. Galactic Center Pulsar Population 

We show here that the pulsar detections within 15 arcmin of Sgr A* imply the existence 
of a Galactic center subpopulation of pulsars and that this population may be quite sizable. 
First, given the small solid-angle coverage of our survey, AQ pa 10~ 4,2 sr, we do not expect to 
find any pulsars from the Galactic disk population. Integrating ei ther a simple disk m odel 



with a uniform density or using the spatial distribution reported by lLorimer et al.l (120061 ). we 
find that only 1 to 2 disk pulsars beamed toward us are expected in the volume between Earth 
and Sgr A*. When we include the survey sensitivity, which depends on the pulsar period, 
and the luminosity and period distributions, the expected number of detected pulsars <C 1. 

To assess the rough parameters of a Galactic center pulsar population, we employed 
a simple modeling approach. We used an ellipsoidal distribution with characteristic ra- 
dius in the plane of the Galaxy Rqc that contains Nqc pulsars. The scale height of the 
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Fig. 2.— Pulse profiles at 1.5, 2, 5, and 9 GHz for J1746-2850I (left), J1746-2850II (middle), 
and J1745-2910 (right). Scattering broadening results in an exponential tail to the pulse 
profiles at 2 and 1.5 GHz. J1746-2850II is scattered most severely of the three pulsars due 
to its high DM. 
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Fig. 3. — Expected dispersion measure vs. distance for sever al lines of sight in th e inner 
Galaxy based on the NE2001 model of Galactic ionized gas (ICordes fc Lazid 120031 ) . The 
inset shows expected DM vs. distance for the lines of sight to the five known pulsars within 
1° of Sgr A*. The expected DM for objects physically near Sgr A* is > 1500 pc cm -3 . 
The vertical part of the curve for a line of sight towards the Galactic center corresponds 
to a den se scattering screen ar ound Sgr A* which makes detection of pulsars in this region 
difficult jLazio fc Corded Il998h . 
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Fig. 4. — Likelihood contours based on the Poisson probability of detecting on average 
(N) pulsars of the simulated population in a survey beam where there are k actual de- 
tections. Contours from left to right and in order of increasing thickness correspond to 
C = (0.1,1,2,4,6,8) x 1(T 6 . 
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population was fixed at Hqq — 0.026 kpc, equal to that of the scattering region in the 
NE2001 model. Through Monte Carlo, we generated pulsars consistent with this spatial 
distribution a n d dra wn from populations with the period and luminosity distributions of 



Lorimer et al.l (120061 ). Detection criteria were applied using the GB T observation parame - 
ters reported earlier and the Parkes survey parameters reported by IJohnston et al.l (120061 ). 



and taking into account the effects of pulse broadening from both residual dispersion smear- 
ing and scattering on the harmonic sum. As a test statistic we used the likelihood function 
C = Y\iP{h,bi, (Ni) ,ki). A term in the product corresponds to the Poisson probability 
of detecting on average (N) pulsars of the simulated population in the z-th survey beam, 



where there are k actual detections, P((N) , k) = (N) k jk\. The population was gener- 
ated 1000 times for each Nqc,Rgc pair. Figure 0] shows contours of C, which indicate that 
iVcc > 2000 and Rqc ^ 0.3 kpc are lower bounds on parameter values. 

Extending our grid up to Nqc = 10 4 and Rqc = 5 kpc, we found no upper bounds on 
the parameters. We expect such bounds from a more comprehensive analysis that includes 
other survey results over a larger region, which we defer to another paper. 



5. Conclusions 



We have discovered three pulsars within 12' of SgrA*. The low (<C 1) number of 
detectable disk pulsars in our survey volume and the high dispersion measures of the new 
pulsars indicate that they are within the dense scattering region surrounding SgrA* and 
part of a neutron star population associated with the Galactic center. Based on a Monte 
Carlo simulation of th is population which incorporates our survey results as well as those of 
Johnston et al.l (120061 ). we conclude that there are Nqc ^ 2000 active pulsars associated with 
the Galactic center. Based on the average lifetime of canonical pulsars, ~ 10 Myr, we obtain 
a rough estimate of the radio pulsar birth rate in the Galactic center of > 2 x 10 -4 yr^ 1 . 
Considering that this is a lower limit, and that not all supernovae produce radio pulsars, 
this result is consistent with the rate of Galactic center su pernovae (~ 10~ 3 yr _1 ) based 



on a survey of co mpact radio sourc es (ILazio fc Corded 120081 ) and observations of soft X-ray 
emitting plasma (IMuno et al.l 12004 ). 
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